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Abstract. We propose that some of the unidentified 
\l Galactic plane (—10° < b < 10°) sources listed in the 



On 
ON 



first EGRET catalog (Fichtel et al. 1994) are associated 
with nearby supernova remnants. The probability that the 
association of strong EGRET sources with supernova rem- 
^^nants is due to chance alignment is found to be ^ 1/2200. 
fyy Three of the most convincing associations are radio-bright, 



nearby shell remnants and a fourth is associated with 



a plerion. We examine possible gamma-ray production 
mechanisms and conclude that pion production and decay 

^ from high-energy protons interacting with the remnant gas 

could produce the gamma rays. The energetics of this pro- 
C^} cess are in accord with theories of cosmic ray production 
ON through supernova shock acceleration. The identification 
of 7r°-decay features in the gamma-ray spectra of super- 



On 



nova remnants would support this interpretation. 
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1. Introduction 

The identified high-energy gamma-ray ( 100 MeV) 
sources observed with the EGRET experiment on the 
Compton Gamma-Ray Observatory have fallen into five 
general categories: solar flares, pulsars, gamma-ray bursts, 
normal galaxies, and active galactic nuclei. In addition, 
the first EGRET source catalog (Fichtel et al. 1994) lists 
37 unidentified Galactic plane sources, of which 10 are 
high-confidence detections (> 6<r) and 27 are marginal 
detections (between 5 and 6<r). We have compared the lo- 
cations of these 37 sources with the locations of the 182 
supernova remnants (SNRs) listed in the catalog of Green 
(1994). For simplicity, the radii of elliptical SNRs were 
taken to be the average of the semi-major and semi-minor 
axes of the remnant. We found that of the 37 unidentified 



high-energy gamma-ray sources, 4 of the high confidence 
sources and 9 of the marginal detections have SNRs whose 
radii overlap the 95% confidence EGRET gamma-ray er- 
ror circle (one of the marginal detections has two possible 
associations; see Table 1). Two of the gamma-ray sources 
have SNRs within 8' of the centroid of the EGRET error 
circle. One of these two sources is at a Galactic longitude 
of 189° where there are few SNRs and the possibility of 
chance coincidence is low. 



2. Statistical Test 

We developed a test to determine if the associations be- 
tween the 37 EGRET sources and the 182 SNRs are chance 
alignments. We calculate the angular distance, as mea- 
sured from the center of the EGRET error circle to the 
center of the remnant, from each gamma-ray source to the 
closest three SNRs, 8±, 62, and #3. Assuming that the near- 
est neighbors are randomly distributed about the EGRET 
source, we then use these distances to define two param- 
eters, ai = {6\/6l) and a 2 = (Q\ - 6\)/{fil - The 
quantity a\ is the ratio of the areas of circles with radii 6\ 
and 62 centered at the location of the EGRET source. The 
quantity 0,2 is the ratio of the areas of annuli with inner 
radius 6\ and outer radii 62 and #3, respectively, centered 
on the EGRET source. If the SNRs are not correlated with 
the EGRET sources, then (ai), the average value of a±, 
should equal 0.5 within statistical errors. If the unidenti- 
fied EGRET sources are, however, associated with SNRs, 
then (ai) should be < 0.5. The average value of 0,2, (0,2)1 
should be 0.5 if our assumption regarding the distribution 
of nearest neighbors is correct. 

Another way of testing the assumption that the distri- 
bution of the nearest-neighbor SNRs about an unidenti- 
fied EGRET source is two-dimensional is to compute the 
average of the sines of the angles between the directions 
from the EGRET source to the nearest and second nearest 
(second and third nearest) SNRs, (sin-0 1 ( 2 ))- If the SNR 
distribution is two-dimensional, (sin-0) = 2/w = 0.637. If 
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Table 1. Unidentified EGRET Sources With Possible SNR Associations 



EGRET Source 


SNR 


0i (T 


n ("I 6 

max \ J 


$1 / D max 


Type c 


Radio Flux (Jy) 


GRO J0542+26 


G 180.0-1.7 (S147) 


116.6 


248.0 


0.47 


s 


65 


GRO J0617+22 


G 189.1+3.0 (IC 443) 


6.7 


43.5 


0.15 


s 


160 


GRO J1110-60 


G 291.0-0.1 (MSH 11-62) 


7.8 


56.0 


0.14 


F 


16 


GRO J2019+40 


G 78 2+2 1 (-v Cvenil 


27.7 


48.0 


0.58 


s 


340 


GRO J0635+05 


G 205.5+0.5 (Monoceros) 


81.6 


148.0 


0.55 


s 


160 


GRO J0823-46 


G 263.9-3.3 (Vela) 


87.7 


127.5 


0.69 


c 


1750 


GRO J1416-61 


G 312.4-0.4 


11.9 


49.8 


0.24 


s 


44 


GRO J1443-60 


G 316.3+0.0 (MSH 14-57) 


35.4 


39.0 


0.91 


s 


24 


GRO J1758-23 


G 6.4-0.1 (W28) 


26.8 


57.0 


0.47 


c 


310 


GRO J1823-12 


G 18.8+0.3 (Kes 67) 


13.0 


36.8 


0.35 


s 


27 


GRO J1842-02 


G 30.7+1.0 


42.3 


50.5 


0.84 


s 


6 


GRO J1853+01 


G 34.7-0.4 (W44) 


30.7 


42.5 


0.72 


s 


230 


GRO J1904+06 


G 40.5-0.5 


36.7 


63.0 


0.58 


s 


11 




G 41.1-0.3 (3C397) 


46.1 


53.8 


0.86 


s 


22 



a Angular distance from center of EGRET error circle to center of associated remnant. 
b Sum of the EGRET error circle radius plus the radius of the associated remnant. 
c S=Shell, C=Composite, F=Filled 



the distribution is linear, the average should be much less 
than 0.637. We find that for the 10 high confidence sources 
(sin-0 1 ( 2 )) = 0.69(0.47) ± 0.10 where the 1<7 error is given 
by (0.5-4/7r 2 ) 1 / 2 / 1 /^ 7 ! where N is the number of sources. 
For the 27 marginal sources, (sin-0 1 ( 2 )) = 0.53(0.64) ±0.06 
and for the entire sample (sin-0 1 ( 2 )) = 0.58(0.60) ± 0.05. 
Thus our assumption appears valid. 

We find that given the entire sample of 37 unidenti- 
fied high-energy gamma-ray sources, (a±) = 0.41 ± 0.05 
and (a 2 ) = 0.45 ± 0.05, where the errors are statistical, 
i.e. it = l/-\/T2N where N is the number of sources. Thus 
there is ~ 1.8(7 deviation of (a±) from 0.5 while (a 2 ) de- 
viates by only 1<7 from 0.5. This weakly suggests that the 
distribution of SNRs closest to the gamma-ray sources is 
not completely random while supporting the assumption 
that the distribution of second closest remnants is ran- 
dom. When only the 10 high confidence sources are stud- 
ied, however, the values of (a±) and (a 2 ) are 0.21+0.09 and 
0.53 + 0.09, respectively. In this case, there is ~ 3.2(7 devi- 
ation from 0.5 for (a±) while (a 2 ) is within 1<7 of 0.5, sug- 
gesting only a small chance of coincidental alignment with 
the closest remnants. Computer simulations show that the 
probability of having (cti) < 0.21 is < 1/2200. 

When considering only the marginal gamma-ray source 
detections, (cti) = 0.48 + 0.06 and (a 2 ) = 0.43 + 0.06, both 
roughly consistent with 0.5. The lack of statistically sig- 
nificant association between the weak EGRET-source and 
supernova-remnant populations is in part because most 
marginal EGRET detections are located toward the galac- 
tic centre region where there is a larger density of SNRs, 
and the relatively large size of the EGRET error circle re- 
duces the probability of association. We show histograms 



of the values of a\ for the high-confidence, marginal, and 
total set of EGRET sources in Figure 1. Note that all but 
one of the high-confidence sources has a value for a\ < 0.5. 
We would like to note that Esposito et al.(1994) have inde- 
pendently found a statistically significant correlation be- 
tween the unidentified EGRET sources and SNRs using a 
different statistical test. 

3. Properties of the Associated Supernova Rem- 
nants 

We now consider the energetics of the SNRs associ- 
ated with the high confidence gamma-ray detections. Of 
these, three are radio-bright, nearby shell-type remnants 
(Green 1994), namely 7 Cygni (G 78.2+2.1), IC 443 (G 
189.1+3.0), and S147 (G 180.0-1.7), and the fourth, MSH 
11-62 (G 291.0-0.1), is a plerion (Weiler &; Sramek 1988). 
The distance to IC 443 is between 0.7 and 2.0 kpc (Green 
1994), and Petre et al. (1988) have estimated its age at 
2800-3400 years. 7 Cygni has an estimated distance of 1.5 
kpc (Green 1989; Landecker, Roger, &: Higgs 1980). An 
age for 7 Cygni can be estimated assuming that this rem- 
nant is in the Sedov phase using the equation (Kassim, 
Hertz, &; Weiler 1993) 

r^O.l {d6) 25 (^-) years, (1) 

where d is the distance to the remnant in kpc, 6 is the 
angular diameter of the remnant in arcminutes, n is the 
density of the interstellar medium surrounding the rem- 
nant in cm -3 , and €51 is the supernova explosion energy 
in units of 10 51 ergs. The ISM density n a typically has val- 
ues between 0.1 and 1.0 cm -3 . For simplicity we choose 
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2 shell remnants. The age of this remnant along with the 
large 95% confidence error circle (158') and the large size 
of the remnant make chance alignment a distinct possibil- 
ity, although it is located in a region with few other rem- 
nants. MSH 11-62 is a plerion with an age of 3000 years 
and a distance of ~ 3.5 kpc (Roger et al. 1986). Therefore 
its gamma rays could originate from high-energy particles 
emitted by the pulsar which is presumably responsible for 
powering this remnant. 

We have ordered Green's (1994) catalog of SNRs ac- 
cording to radio flux. 7 Cygni and IC 443 are the fifth 
(340 Jy) and eleventh (160 Jy) brightest SNRs at 1 GHz, 
respectively. S147 and MSH 11-62 are significantly dim- 
mer in the radio, and are numbers 25 (65 Jy) and 74 (16 
Jy) on the list. When compared with only the flatter spec- 
trum (spectral index ;$ 0.5) shell remnants, 7 Cygni and 
IC 443 are the second and seventh brightest. The radio 
bright SNRs that are not associated with EGRET sources 
tend to be either farther away than 7 Cygni and IC 443, 
such as Cas A and W28 at 2.8 and 3.5-4.0 kpc (Green 
1994), respectively, or are much older, such as the Lu- 
pus Loop at ~ 50, 000 years (Leahy, Nousek, &: Hamilton 
1991). 



Fig. 1. A histogram of the number of EGRET unidentified 
sources with different values of the parameter ai = 
where Q\ and 02 are the angular distances to the closest and 
second closest supernova remnants, respectively. Note the large 
fraction of high confidence EGRET sources which have values 
of ai < 0.1, indicating likely associations. 



n /e 51 = 0.5. Thus r k, 0.07 (d 0) 2 5 years. For 7 Cygni, 
t 5400 years given its 60' diameter (Green 1994). For 
comparison, Equation (1) gives an age of 13,000 years for 
the Vela supernova remnant which is in good agreement 
with the 11,000 year spin-down age for the Vela pulsar. On 
the other hand, S147 is a nearby, old shell-type remnant. 
Sauvageot, Ballet, &: Rothenflug (1990) give its distance 
as 1.1 — 1.7 kpc. Given its 180' diameter, Equation (1) 
yields an age of ~70,000 years, much older than the other 



4. Possible Gamma- Ray Production Mechanisms 

Gamma-ray emission from MSH 11-62 can be readily un- 
derstood as either pulsed emission from an associated 
radio-quiet pulsar with unknown period or as unpulsed 
emission from a pulsar-energized nebula (plerion) such as 
the Crab Nebula. Detection of a point source in X-rays by 
ROSAT or ASCA would give evidence for a pulsar source 
for the emission and could permit a search for the pul- 
sar's period. Gamma rays from shell-type remnants such 
as 7 Cygni and IC 443 (and possibly S147) may involve 
other production mechanisms, although searches for radio- 
quiet pulsars associated with these sources should also be 
pursued. Let us therefore consider the possible origin of 
the gamma rays from the shell-type SNRs. It is not likely 
that the high-energy radiation originates from Compton 
upscattering of a photon field by the synchrotron-emitting 
electrons. For example, IC 443 at an assumed distance of 
1.5 kpc, emits k, 10 33 ergs s _1 in the radio (Green 1994; 
Mufson et al. 1986), « 10 37 ergs s _1 in the IR (Mufson 
et al. 1986), and 10 35 ergs s -1 in X-rays (Petre et al. 
1988; Mufson et al. 1986). If our proposed association is 
correct, IC 443 emits « 5 x 10 34 ergs s _1 in > 100 MeV 
gamma rays (Fichtel et al. 1994). The relationship between 
the gamma-ray luminosity L 1 and the radio luminosity L r 
from Compton upscattering and synchrotron radiation, re- 
spectively, can be shown to be 



'->ph 



h. ~ 50 

L r 47rr 2 (5 2 /87r)c' 



(2) 



where L p h is the luminosity of one of the lower-energy 
photon fields, B is the mean magnetic field strength in 
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the remnant shell [typically of order 10 -5 Gauss (Lozin- 
skaya 1992)], and r is the remnant radius (the rem- 
nant thickness does not enter). Equation (2) implies that 
5(10" 5 Gauss) « 0.01(£ pft /10 37 ergs s' 1 ) 1 / 2 /(r/10pc). 
Thus, unless the remnant magnetic field is extremely 
weak, there are not enough synchrotron-emitting electrons 
to produce the observed gamma rays by upscattering the 
known radiation fields or, for that matter, the microwave 
background or starlight radiation fields. 

The relativistic electrons which produce the observed 
radio synchrotron radiation will also produce gamma- 
rays through nonthermal bremsstrahlung. The nonther- 
mal electron spectrum in IC 443 can be estimated by 
comparing the derived synchrotron photon luminosity 
with its observed radio luminosity. The synchrotron or 
bremsstrahlung luminosity can be written as 

L = m e c 2 N e (j e )jdj e , (3) 

where, for synchrotron radiation, j = j syn — 
(4/3)c<7 TTe 2 (5 2 /87rm e c 2 ) and N e (j e ) = N°y-i is the elec- 
tron density per unit j e . In order to produce radio photons 
with frequencies near 10 GHz, To > 1.9 x 10 4 5I°- 5 , where 
5_5 is the magnetic field in units of 10 -5 Gauss. The 
spectral index of the electron population, q, must Ri 1.72 
in order to produce the 0.36 radio spectral index observed 
for IC 443 (Green 1994). Setting Equation (3) equal to 
10 33 £ 3 3 ergs s" 1 , we find N° « 4 x 10 52 5Z^ 36 £ 3 3- For 
bremsstrahlung radiation, j = jbrem — (l/3) c<7 T£*/ra/j7 e , 
and tih = lOrai cm -3 is the mean number density of hy- 
drogen atoms in the remnant. The derived bremsstrahlung 
luminosity Ri 9 x 10 32 niBZl' 5 L33 ergs s _1 , which im- 
plies 5(10 -5 Gauss) k, 0.07(7Xi£33) - 67 given the observed 
gamma-ray luminosity of IC 443. 

It was suggested long ago that if SNRs are sites of 
cosmic-ray acceleration, then 7r°-decay gamma-rays will 
be produced as a result of collisions of the cosmic rays with 
the interstellar medium (Pinkau 1970). The luminosity in 
secondary 7r°s is 

/•OO 

Lt° — cn H I N p (jp)o-pp(jp)(E(jp))dj p , (4) 

" 7p,thresh~2 

where N p (j p ) = N°j~ 2 ' 2 is the number of protons per unit 
j p suggested by shock acceleration models (Gaisser 1990), 
a pp 3 x 10 _26 cm -2 (Dermer 1986) is the proton-proton 
interaction cross section, and (E(j p )) = l/3(7 p m p c 2 /2) is 
the average amount of energy per proton-proton interac- 
tion converted to 7r°-decay gamma rays. By setting Equa- 
tion (4) equal to the gamma-ray luminosity of the EGRET 
source possibly associated with IC 443, k, 5x 10 34 ergs s _1 , 
we find N° k, 5 x 10 51 n^ 1 protons. Thus the amount of 
energy in nonthermal protons required to produce the ob- 
served gamma-ray emission is 4 x 10 49 n^ 1 ergs. This 
can be compared with E cr , the amount of energy per su- 
pernova that is converted to nonthermal protons in order 



to produce the observed cosmic ray energy density ty cr , 
which we write as E cr = ^ C rVg a ir sn /r ga i. Here * cr = 
10" 12 *_i2 ergs cm" 3 (Gaisser 1990), V gal = 10 67 V 6T cm 3 
is the volume of the galaxy, r sn = 10 9 t sn ^ seconds is the 
time between supernovae in the galaxy, r ga \ = 3 x 10 14 t? 
seconds is the proton confinement time in our galaxy, 
and t-j is the confinement time in units of 10 7 years. 
Thus E cr Ri 3 x 10 49 ^ -i2Vert S7lj g/tr ergs, similar to that 
needed to produce the observed gamma-ray flux. Aharo- 
nian, Drury, and Volk (1994) have also found that gamma 
rays produced by the decay of 7r°s could be detectable with 
current instruments, particularly if the remnant is adja- 
cent to a molecular cloud as is evidently the case for 7 
Cygni and IC 443 (Pollock 1985; Asaoka &: Aschenbach 
1994) Observations of pion-decay features at gamma-ray 
energies from SNRs with the telescopes on the Compton 
Observatory would support the w° origin of the gamma-ray 
emission and be in accord with the scenario that cosmic 
rays are accelerated by supernova shocks. 
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Figure 1 



